Observation of unconventional band topology in a superconducting doped topological 
insulator, Cu x Bi 2 Se 3 : Topological or non-Abelian superconductor? 
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The Cu-doped topological insulator Bi 2 Se3 has recently been found to undergo a superconducting 
transition upon cooling, raising the possibilities that it is the first known "topological superconduc- 
tor" or realizes a novel non-Abelian superconducting state. Its true nature depends critically on the 
bulk and surface state band topology. We present the first photoemission spectroscopy results where 
by examining the band topology at many different copper doping values we discover that the topo- 
logically protected spin-helical surface states remain well protected and separate from bulk Dirac 
bands at the Fermi level where Copper pairing occurs in the optimally doped topological insulator. 
The addition of copper is found to result in nonlinear electron doping and strong renormalization 
of the topological surface states. These highly unusual observations strongly suggest that supercon- 
ductivity on the topological surface of CuccBi2Se3 cannot be of any conventional type in account of 
the general topological theory. Characteristics of the three dimensional bulk Dirac band structure 
are reported for the first time with respect to the superconducting doping state and topological 
invariant properties which should help formulate a specific theory for this novel superconductor. 

PACS numbers: 



Topological insulators embody a new state of matter 
characterized by topological invariants of the band struc- 
ture rather than spontaneously broken symmetry, and 
feature massless Dirac-like conduction states on their sur- 
faces [l], Q, H, 0] • Bismuth selenide in particular has been 
found to be an ideal "hydrogen atom" topological insu- 
lator, realizing the simplest known case of tojjologically 
nontrivial band structure 

H % SB G|. It has 

been proposed that inducing a superconducting gap in 
the surface states of topological insulators will lead to 
fault tolerant non-Abelian surface physics with potential 
application in spintronics and quantum computing [T3 |. 
In this Letter, we use angle resolved photoemission spec- 
troscopy (ARPES) to examine the effect of copper dop- 
ing on the electron dynamics of Cua3i2Se3, which has 
been shown to lead to a bulk superconducting transition 
at x>0.1 (max Tc=3.8° K) [16[. Band structure in the 
normal state superconductor is found to preserve a min- 
imalist scenario for a topological metal, consisting of a 
single bulk band with nearly isotropic massive Dirac-like 
dispersion and a well defined, topologically protected sur- 
face Dirac cone. As a result of this band structure, we 
find that superconductivity at the sample surface can- 
not be conventional. Copper doping is systematically ob- 
served to have complex effects on the system, adding a 
small number of electrons and causing a strong renormal- 
ization of the surface bands in such a way as to preserve 
the topological character. The doped compound exhibits 
hexagonal dispersion anisotropy, with important impli- 
cations for low energy interactions and potential device 
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FIG. 1: Superconductivity in a doped topological in- 
sulator Bi2Se3: (a) Resistivity and magnetic susceptibility 
measurements from Ref. [16|] show a superconducting transi- 
tion temperature of 3.8°K at optimal copper doping (x=0.12). 
(b) The hexagonal surface state Brillouin zone of Cu ;E Bi2Se3 
is drawn in red above a diagram of the three dimensional bulk 
Brillouin zone, (c) A LEED image taken at 200eV shows a 
well ordered surface with no signs of superstructure. 



development [13[. 

Undoped Bi 2 Se3 is a topological insulator with a large 
band gap (>300 meV) [5j, and belongs to a class of mate- 
rials M2X3 (M=Bi,Sb; X=S,Se,Te) that includes at least 
two other topologically nontrivial materials, Bi2Te3 and 
Sb2Te3, with smaller band gaps and more complicated 
band structure [1, 0] . These materials share a rhombo- 
hedral crystal structure, with a five atom unit cell ar- 
ranged in quintuple layers, and have been investigated 
extensively in connection to thermoelectric applications 
[lij . Unlike band structure in topologically trivial ma- 
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terials, which is commonly more parabolic ("classical" 

2 

E=^g), it is most natural for the bulk conduction bands 
of topological insulators to realize Dirac-like dispersion 
following an analogue of Einstein's equation for the en- 
ergy of a relativistic particle (E 2 =M 2 v^+p 2 v^) with a 
critical velocity "vc" analogous to the speed of light. 
The addition of electrons from copper doping allows us 
to quantitatively evaluate the bulk Dirac-like band char- 
acter in the relativistic regime. 

Angle resolved photoemission spectroscopy (ARPES) 
measurements were performed at the Advanced Light 
Source beamline 10.0.1 using 35.5-48 eV photons and 
Stanford Synchrotron Radiation Laboratory (7-22eV 
photons) with better than 15 meV energy resolution and 
overall angular resolution better than 1% of the Brillouin 
zone (BZ). Samples were cleaved and measured at 15°K, 
in a vacuum maintained below 8xl0 -11 Torr. Momen- 
tum along the z-axis is determined using an inner poten- 
tial of 9.5 eV, consistent with previous ARPES investi- 
gations of undoped Bi2Se3 [5[. 

Surface and bulk state band calculations were per- 
formed for comparison with the experimental data, using 
the LAPW method implemented in the WIEN2K pack- 
age [Hj. Details of the calculation are identical to those 
described in Ref. |5[. 

To better understand the highly nonlinear doping ef- 
fect of copper, we present data on several copper-added 
Cu^Bi2Se3 crystals (x=0, 0.01, 0.05, 0.12) and copper 
substituted Cuo.1Bi1.9Se3, grown as described in Ref. 
[l6j |. Copper atoms intercalated between the van der 
Waals bonded selenium planes are thought to be sin- 
gle electron donors, while substitutional defects in which 
copper replaces bismuth in the lattice (Cu#i) are ex- 
pected to each contribute two holes to the system (re- 
introducing x=0.12 copper doping for optimal supercon- 
ductivity has been found to shift the z-axis lattice pa- 
rameter by only 1.5% while leaving the in-plane lattice 
parameters and long-range crystalline order intact [16j . 

The qualitative effect of copper addition is an enlarge- 
ment of the Fermi surface from electron doping and, sur- 
prisingly, a strong renormalization of the surface state. 
The surface state dispersion is flattened by 30% when 
copper doping of x=0.05 is added to the stoichiometric 
compound, and the Fermi surface becomes hexagonally 
anisotropic in a way consistent with numerical predic- 
tions in the local density approximation (LDA). Tunabil- 
ity of the surface state anisotropy is important for con- 
trol of unconventional ordered states that may appear 
uniquely in topologically ordered materials [13[ . 

The bulk and surface-derived band structure of the 
normal state optimally doped superconductor (x=0.12) 
are explored in Fig-3. The bottom of the conduction 
band for electrons in the bulk is found at the three dimen- 
sional T-point (k z =k x =k y =0) , and can be seen inside the 
upper surface state (SS) cone. Due to the renormaliza- 



tion effect and less-than-linear electron doping, the occu- 
pied surface state band structure does not intersect with 
the bulk conduction band structure at any point in mo- 
mentum space, making the surface state "well defined" 
and topologically protected. The bulk conduction band is 
only clearly visible when low photon energies (h^<20eV) 
are used to increase bulk penetration, probably due to a 
screening effect related to excess negative charge carriers 
in the surface state, (see analysis in Fig-2(d)) The gap 
between bulk valence and conduction bands appears to 
be unchanged upon copper doping. 

Using 9.75eV photons to view the T-M and T-K direc- 
tions shows Fermi momenta of 0.110±3 A -1 and 0.106±3 
A -1 respectively. Varying incident energy to observe dis- 
persion along the z axis (T-Z direction) reveals a Fermi 
momentum of 0.12 A -1 , suggesting that the bulk elec- 
tron kinetics are three dimensionally isotropic. Carefully 
tracing the band (Fig-3 (d-e)) yields a Fermi velocity of 
3.5 eVxA along T-M and 4.1 eVxA along T-K, esti- 
mated within 50 meV of the Fermi level. Assuming that 
the bulk conduction band forms a Fermi sea for super- 
conductivity, the superconducting correlation length can 
be estimated based on the average Fermi velocity and 
superconducting critical temperature to be about 2000A 
(£o~ 0.2 x hv F /K B T c = 0.2 x 3.8eV A/(K B x 3.8°K) = 
2000 A), large enough for phase fluctuations to be ne- 
glected in the neighborhood of Tc*. This is a typical 
coherence length for conventional superconducting ma- 
terials, and much greater than that seen in other uncon- 
ventional superconductors, such as strongly correlated 
cuprates (£ ~ 100-200A) or cobaltates (f ~ 200A) 0. 

Along both the T-M and T-K directions, dispersion of 
the bulk conduction band appears to approach a limit- 
ing velocity at large momentum, rather than following a 
parabolic arc. This behavior, like the "massless" linear 
dispersion of the surface state bands, is mathematically 
analogous to the energy-momentum relationship of a rel- 
ativistic particle approaching the speed of light. Classi- 
cal (paraboloic, m=0.155 m e ) and relativistic Dirac-like 
(m=0.155 m e , v c =6 eV-A) energy dispersions are plotted 
as fits for the bulk band in Fig-3 (e). We find that a rel- 
ativistic fit can more accurately reproduce the observed 
dispersion, however neither fit can account for a slight 
bend in the dispersion centered near 90 meV, which may 
suggest strong electron-phonon interactions in the sys- 
tem consistent with phonon-mediated superconductivity. 
The simplicity of band structure in the doped Bi2Se3 
system clearly sets it apart from other known topological 
insulators (e.g. Bi2Te3) that have far greater deviation 
from pure Dirac-like kinetics @, 0, EH, El . 

The weakness of copper doping is critical to preserv- 
ing topological order in the superconducting state, by 
keeping the Fermi level beneath the point of intersection 
between bulk and surface conduction bands, (near the 
top of the diagram in Fig-4(a,right)) The number of con- 
ducting charge carriers at the surface and in the bulk are 
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FIG. 2: Nonlinear doping and band topology: (a) Symmetrized Fermi surfaces are displayed together with (right) an 
LDA prediction based on the Dirac point energy at x=0.05 copper doping. For clarity the middle panel is measured at a 
photon energy that suppresses bulk resonance (48eV). (b) ARPES measurements are shown at high photon energy (E>20 eV) 
for non-superconducting Cu x Bi2Se3. Particle velocities in the surface state upper Dirac cone are reduced by 30% after the 
addition of x=0.05 copper doping, (c) When the Bi-Se plane spacing at the surface of a 12-layer slab is increased by 0.2 A (1.64 
to 1.84 A), dispersion in the upper Dirac cone increases by 16%, similar to the amount of renormalization we have observed in 
the superconducting compound. The surface state Dirac point is set at zero energy for each calculation, (d) The total number 
of charge carriers in the bulk and at the surface is calculated from the Luttinger count ( ^ , x2 for the doubly degenerate 
bulk band). 




FIG. 3: Bulk Dirac states vs Surface band topology and superconductivity: (a-b) Momentum dependence of the 
bulk and surface conduction bands in superconducting Cuo.i2Bi2Se3 is measured through the 3D Brillouin zone center with low 
energy (9.75eV) photons for enhanced bulk sensitivity, (c) Small panels show that the bulk and surface bands remain separate 
at intermediate k z values, (d) Dispersion is traced on a close-up image of the conduction bands, (e) Energy-momentum 
dispersion of the bulk electrons is compared with Dirac-like (v c =6 eV-A) and classical (parabolic) fits with an effective mass 
of 0.155m e . Neither fit can perfectly account for a slight bend in dispersion near 90 meV binding energy. An inset shows the 
difference between bulk band energy and the Dirac-like fit curve. (f,blue) The bulk band dispersion approximated from ARPES 
data is plotted at k z =0, surrounded by (green) a massless "light-like" limit, illustrating that bulk electron kinetics are in the 
"relativistic" Dirac regime. 



estimated using the Luttinger count in Fig-2(d), by divid- 
ing the surface and bulk Fermi surface areas by the total 
area of the Brillouin zone. Although the doping increases 
monotonically as copper is added, the carrier density in 
the bulk is only ~l/30th of what would be expected if all 
copper were intercalated between paired selenium layers. 
It is therefor likely that nearly one third of the added 
copper enters the sample through substitutional defects 



with bismuth, adding holes that counterbalance most of 
the electron doping from intercalation. 

Attempting to force the creation of Cu#i defects by 
adding less bismuth results in very weak hole doping for 
Cuo.1Bi1.9Se3. In this case, the binding energy of the 
Dirac point decreases by more than 100 meV relative 
to the undoped compound, raising the bulk conduction 
band entirely above the Fermi level so that the material 
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Band structure based on experimental data 




FIG. 4: Symmetry breaking and the topological sur- 
face: (a) Topologically protected surface states cross the 
Fermi level before merging with the bulk valence and conduc- 
tion bands in a topological material, (b) If superconducting 
parity is even, the surface states will be gapped due to their 
proximity to the bulk superconducting wavefunction, allowing 
non-Abelian Majoran Fermion surface vortices, (c) If parity 
is odd the material will be a topological superconductor, and 
new states will appear below Tc to span the bulk supercon- 
ducting gap. (d,e) States within 5meV of the Fermi level are 
shaded for the k 2 =0 plane of (Case 1) a fully gapped doped 
topological insulator and (Case 2) one example of a topolog- 
ical superconductor. 

is a traditional topological insulator. An explanation for 
why hole doping is not stronger in this case may be that 
defects in which selenium fills bismuth vacancies will re- 
sult in the addition of electrons. 

The surface state renormalization that we have ob- 
served is also instrumental in separating the bulk and 
surface band structure. Renormalization is likely to re- 
sult from a combination of increased carrier density on 
the surface relative to the bulk, evident from Fig-2(d), 
and relaxation of interlayer bond lengths near the cleaved 
surface. Changes to interatomic bonding lengths at the 
surface of strongly spin-orbit coupled systems can have 
a significant effect on surface state dispersion [20[, and 
could be effected by the surface carrier density and pres- 
ence of copper. We find that increasing the distance be- 
tween the outermost bismuth and selenium layers by 0.2 
A can account for the difference in dispersion between 
undoped and superconducting doped crystals (Fig-2(c)). 

A cartoon of band structure in a topological metal (in- 
sulator) is overlaid on ARPES data in Fig-4(a), illustrat- 
ing that the gapped bulk bands are connected by a singly 
degenerate gapless surface state. When a bulk super- 
conducting transition is introduced, there are two likely 
scenarios. If the parity eigenvalue of the bulk supercon- 
ducting state is even (Fig-4(b,d)), electrons in the surface 
state will participate in bulk superconductivity through 
the proximity effect, resulting in the appearance of non- 
Abelian vortex states of great interest for quantum com- 
puting [12[. If parity is odd, the system will be a topolog- 
ical superconductor [2l[ , with new gapless states appear- 



ing beneath the Tc. ARPES measurements cannot re- 
solve the bulk superconducting gap, which is expected to 
be only -0.6 meV from BCS theory (3.5xK B T c /2=0.6 
meV). 

Although our results show that Bi2Se3 is a nearly ideal 
minimalist topological insulating system, many of the 
most interesting physical properties emerge due to its 
deviations from the simplest case of perfectly isotropic 
Dirac-like electron kinetics. Recent theoretical explo- 
rations have suggested that magnetic perturbations will 
not readily open a mass gap in the surface state of a topo- 
logical insulator with perfectly linear dispersion [22|, , 
but can do so when the bands have some upward con- 
cavity as in Bi2Se3 [2J]. The hexagonal Fermi surface 
anisotropy demonstrated upon copper doping adds an 
out of plane component to spin polarization, and is a 
likely precursor to novel two dimensional ordered states 
unique to the surfaces of materials with topological or- 
der [l3[. Furthermore, our observation of large surface 
charge contained within the topologically protected sur- 
face state and a dramatic band renormalization effect 
have important implications for manipulation of topo- 
logical insulator systems in future experiments and for 
technological applications. 

In summary, we have investigated the surface and bulk 
band topology of the superconducting doped topological 
insulator Cu :E Bi2Se3. We discover that the copper-doped 
superconducting state unexpectedly preserves the band 
structure and topological order of undoped bismuth se- 
lenide, and that superconductivity at the sample surface 
cannot be conventional. The bulk conduction band con- 
forms approximately to a massive Dirac dispersion, and 
has a slight bend near 90 meV that may be indicative of 
electron-boson coupling. We also observe strong doping 
dependent renormalization of surface states and hexag- 
onal anisotropy that preserve the topological order via 
a 7r Berry's phase invariant and are of interest with re- 
spect to spin transport properties and potential device 
development. 

We thank helpful conversation with L. Fu, B. A. 
Bernevig and P.W. Anderson. 
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FIG. 5: Superconductivity in a doped topological in- 
sulator Bi2Se3: (a) Resistivity and magnetic susceptibility 
measurements from Ref. [16|] show a superconducting transi- 
tion temperature of 3.8°K at optimal copper doping (x=0.12). 
(b) The hexagonal surface state Brillouin zone of CuxE^Ses 
is drawn in red above a diagram of the three dimensional bulk 
Brillouin zone, (c) A LEED image taken at 200eV shows a 
well ordered surface with no signs of superstructure. 
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FIG. 6: Nonlinear doping and band topology: (a) Sym- 
metrized Fermi surfaces are displayed together with (right) an 
LDA prediction based on the Dirac point energy at x=0.05 
copper doping. For clarity the middle panel is measured 
at a photon energy that suppresses bulk resonance (48eV). 
(b) ARPES measurements are shown at high photon energy 
(E>20 eV) for non-superconducting Cu^E^Ses. Particle ve- 
locities in the surface state upper Dirac cone are reduced by 
30% after the addition of x=0.05 copper doping, (c) When 
the Bi-Se plane spacing at the surface of a 12- layer slab is 
increased by 0.2 A (1.64 to 1.84 A), dispersion in the upper 
Dirac cone increases by 16%, similar to the amount of renor- 
malization we have observed in the superconducting com- 
pound. The surface state Dirac point is set at zero energy 
for each calculation, (d) The total number of charge carriers 
in the bulk and at the surface is calculated from the Luttinger 
count ( area ' x2 ^ OY the doubly degenerate bulk band). 
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FIG. 7: Bulk Dirac structure vs Surface band topol- 
ogy and superconductivity: (a-b) Momentum dependence 
of the bulk and surface conduction bands in superconduct- 
ing Cuo.i2Bi2Se3 is measured through the 3D Brillouin zone 
center with low energy (9.75eV) photons for enhanced bulk 
sensitivity, (c) Small panels show that the bulk and surface 
bands remain separate at intermediate k z values, (d) Disper- 
sion is traced on a close-up image of the conduction bands, 
(e) Energy- momentum dispersion of the bulk electrons is com- 
pared with Dirac-like (v c =6 eV-A) and classical (parabolic) 
fits with an effective mass of 0.155m e . Neither fit can per- 
fectly account for a slight bend in dispersion near 90 meV 
binding energy. An inset shows the difference between bulk 
band energy and the Dirac-like fit curve. (f,blue) The bulk 
band dispersion approximated from ARPES data is plotted 
at k 2 =0, surrounded by (green) a massless "light-like" limit, 
illustrating that bulk electron kinetics are in the "relativistic" 
Dirac regime. 
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FIG. 8: Symmetry breaking and the topological sur- 
face: (a) Topologically protected surface states cross the 
Fermi level before merging with the bulk valence and conduc- 
tion bands in a topological material, (b) If superconducting 
parity is even, the surface states will be gapped due to their 
proximity to the bulk superconducting wavefunction, allow- 
ing non-Abelian Majorana Fermion surface vortices, (c) If 
parity is odd the material will be a topological superconduc- 
tor, and new states will appear below Tc to span the bulk 
superconducting gap. (d,e) States within 5meV of the Fermi 
level are shaded for the k 2 =0 plane of (Case 1) a fully gapped 
doped topological insulator and (Case 2) one example of a 
topological superconductor. 



